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TECHNICAL REPORT 



Mostafa Ghandehari 
Mathematics Department 



An Optimal Control Formulation of the 
B 1 aschke - Lebesgue Theorem 



Abstract 



The Blaschke-Lebesgue theorem states that of all plane sets of 
given constant width the Reuleaux triangle has least area. The area 
to be minimized is a functional involving the support function and 
the radius of curvature of the set. The support function satisfies a 
second order ordinary differential equation where the radius of 
curvature is the control parameter. The radius of curvature of a 
plane set of constant width is non-negative and bounded above. Thus 
we can formulate and analyze the Blaschke-Lebesgue theorem as an 
optimal control problem. 
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I nt roduct i on 



The width of a closed convex curve in a given direction is the 
distance between two parallel supporting lines perpendicular to that 
direction. A set of constant width b has the same width in all 
directions. Besides circle, the best known closed convex curve of 
constant width b is the Reuleaux triangle of width b, i.e. a set in 
R 2 whose boundary consists of three congruent circular arcs of radius 
b. See Figure 1. The Blaschke-Lebesgue theorem states that the 
Reuleaux triangle has the least area of all plane convex sets of the 

/_ i ^ \ j^2 

same constant width b. The minimum area is ^ • This theorem 

was first proved independently by Blaschke [2] , and Lebesgue [17] . 
Besicovitch [1] , Chaker ian [9] . and Eggleston [10; 11] contain a proof 
of the Blaschke-Lebesgue theorem. 

Optimal control theory can be applied to geometric extremum 
problems for plane curves as follows: The functional for which 

extrema are examined are geometric invariants such as area or 
perimeter. The system of ordinary differential equations for the 
control theory formulation is derived from the Frenet-Ser ret 
formulas, and the control parameter is curvature. 

Klotzler [16] has used optimal control theory to study n- 
orbiforms. These are convex planar domains which can be rotated 
inside a regular n-gon under tangential contact on all sides. Plane 
sets of constant width can be rotated inside a square with tangential 
contact on all sides. Our approach to plane sets of constant width 
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is different than that in Klotzler in the choice of the functional to 



be minimi zed . 

In the following we discuss preliminary definitions related to 
sets of constant width and include necessary background from the 
theory of optimal control. We then formulate and analyze the 
Blaschke-Lebesgue theorem as an optimal control problem. 

Prel iminaries 

By a convex body in R n we mean a compact convex subset of R n with 
nonempty interior. For each direction u E S"* 1 , where S n_1 is the unit 
sphere centered at the origin in R n , we let h(K,u) denote the support 
function of the convex body K evaluated at u. Thus, 

(1) h(K,u) = s u p { u • x : x E K } , 

which may be interpreted as the distance from the origin to the 
supporting hyperplane of K having outward- po i nt i ng normal u. The 
width of K in direction u is given by 

(2) W ( K , u ) = h ( K , u ) + h ( K , - u ) . 

A convex body K is said to have constant width b if and only if 
W ( K , u ) = b, for all u € S"' 1 . 

For a plane convex body K we shall use the notation h(K,0) = 

h(K,u), where u = (cos 0,sin 0) . In this case equation (2) can be 
written as 

(3) W(I<,0) = h(K,0) + h(K ,0 +tt) . 
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A result that we shall find useful is the formula of Cauchy for the 
Euclidean length of K, namely 

i ( 2ir 

(4) L (K) = A W(K .0) dO . 

J 0 

From (4) we can obtain Barbier J s theorem which states that all plane 
sets of constant width b have the same perimeter 7rb . An elementary 
proof of Barbier’s theorem is given in Honsberger [15] . 

The area of K is denoted by A(K) and is given by 

i ( 2n 

(5) A (K) = A h (K . #) p(K ,0 ) dO 

J 0 

where p(K,0) denotes the radius of curvature of K at the point with 
outward normal u = (cos 0, sin 6) . As a consequence of the Frenet- 
Serret formulas we can obtain the following second order ordinary 
differential equation involving radius of curvature p(K,0) and the 

support function: 

(6) h ( K « 0 ) + h(K,0) = p(K,0) , 

where dot denotes differentiation with respect to 6 . Detailed 
discussion of (4) , (5) and (6) can be found in Flanders [12] , or the 

monograph on convexity by Bonnesen-Fenche 1 [5] . 

For a plane convex body of constant width b we can use (3) to 
obtain 

(7) h(K,0) + h(K,0+? r) = b . 

Substituting 6 by 0+7 r in (6) , we obtain 

(8) h(K,0+7r) + h(K,0+7r) = p(K,6+n) 

Using (7) we have 
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(9) 



h (K , 0) + h (K , 0+ir) = 0 . 

Adding both sides of (6) and (8) and using (9) we conclude that for a 
plane convex body of constant width b the radius of curvature 
sat i sf i es 

(10) p(K,0) + p(K,0+*) = b . 

Since the radius of curvature p(K,0) is nonnegative for a plane 
convex curve, we use (10) to obtain (11) for a plane convex body of 
constant width b. 

(11) 0 < p(K,0) < b . 

The idea of using optimal control is simply to minimize the area 
A(I<) given by (5) subject to differential equation (6) and conditions 
(7) , (10) and (11) . 

In R 3 , the generalization of Barbier 5 s theorem is the fact that 

the total mean curvature M(K) for a set of constant width b is 

constant. That is 

(12) M (I<) = 2?rb . 

We also have the following remarkable relationship for a set of 
constant width b: 

(13) 2V(I<) = bS ( K ) - ^ b 3 , 

where V(K) denotes the volume of K and S(K) the surface area. 
Chakerian [8] contains derivation of (12) and (13). The problem of 
minimizing the volume of a convex body of constant width in R is 
apparently unsolved. Using (13), minimizing the volume is equivalent 
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to minimizing the surface area. The analogues of formulas (5) and 
( 6 ) for volume and the support function of a convex body in R 3 are 
given in Blaschke [3] . The support function H of a convex body in R 3 
sat i sf i es 

(14) AH 4 - 2H = R x + R 2 , 

where R x and R 2 are principal radii of curvature and A is essentially 
the Laplacian restricted to the unit sphere. See Bonnesen - Fenche 1 

[5] . 

Chakerian [9] gives "the following lower bound for n-d i mens i onal 
volume V(K) of a convex body K of constant width 1 i n R n : 



(15) 



V(K) 




k_) 

2 k+ 2 )' 



n > 3 . 



where w n is the volume of the unit ball in R n and A = ^= 7 ^ — . The 

Z 7T 

volume of n -d i mens i onal unit ball in terms of the Gamma function is 
given by 



(16) 



n 

_2 



r (l + §) 



Chakerian and Groemer [7] give an interesting survey article on 
sets of constant width. Bonnesen and Fenchel [5] , Eggleston [11] « 
and Yaglom and Boltyanskii [20] have good treatment of required 
background material from the theory of convex sets. These books also 
contain properties of sets of constant width. 
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Optimal Control 



In the following, we sketch the theory of optimal control and 
give Pontryagin ; s maximum principle. Boltyanskii [4] , Leitmann [18] , 
and Pontryagin [19] are some interesting books on the theory of 
optimal control. Gelfand and Fomin [13] also contains an 
introduction to optimal control. Hermann [14] , and Brockett [6] 
contain differential geometric treatments of calculus of variations 
and control theory. 

We will consider control processes that can be described by a 
system of ordinary differential equations 
dx: 

(1 O cTt = ( X 1 ’ * * * ' x n ’ u l ’ • • • > u r) ? i — 1 5 2 , . . . , n 



where x x , . . . ,x n are space coordinates which characterize the process 
and u 1 ,u 2 , . . . ,u r are control parameters which determine the process. 

In order to determine the process in a given time interval 
[tQjtj] , it is sufficient to give the control parameters as “functions 
o“f time on this interval, that is 

(18) Uj = Uj(t), j = . 

Assuming that the problem is well posed, tor a given initial 
state 

(19) Xj(t 0 ) = Xj° , 



the so 1 ut ion 
Consider 

( 20 ) 



of (17) is uniquely determined, 
the functional 
r * 1 

J = f 0 (x p . . . ,x n ,Up . . . ,u r ) dt . 

J t Q 
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For each control (18) on [to,^] , the process is determined and the 
functional J assumes a certain value. Assuming that there is a 
control (18) which transfers the object from a given initial state 
(19) to a final state 

(21) XjCti) = Xj 1 , 

the object is to find 



(22) Uj(t) , j = 1 r 

which transfers the object from (19) to (21) in such a way that the 
functional (20) has a minimum. In general there are restrictions on 
control parameters Uj . Thus we shall assume the vector u belongs to 
a region U in 2-dimensional Euclidean space called the control 
region, that is 

(23) u G U . 



An admissible control is an arbitrary piecewise continuous 
control in the control region U. 

We now state Pontryagin ’s maximum principle. As a good start for 

understanding, one may study the shortest time problem for a phase 

point moving from an initial point x° to the origin in accordance 

with ^ = x 0 < — = u, lul <1. This problem is studied in many 
dt ^ dt ii_ 

books including Boltyanskii [4] . 



M ax imum Principle 



Let u(t), t Q < t < t x , be 
corresponding trajectory x(t) 



an admissible control such that 
which begins at the point x° at 



the 

the 
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In order that, u(t) 



time t Q passes at the time t, through a point x 1 . 
and x(t) be optimal it is necessary that there exists a nonzero 
continuous vector function ^(t) = (^o(^) ? • • • jV’nC't)) corresponding 
u(t) and x(t) satisfying 



(24) 



d v , i _ m 

dt <9x j 



i — 0,l,2,...,n, 



where 



(25) 

such that 



H(^,x,u) = £ V’afa > 
a=0 



(a) for every t, t Q < t x , the function H(^(t) ,x(t) ,u) of the 
variable u E U attains its maximum at the ponit u = u(t) : 

(26) H(^(t) ,x(t) ,u(t)) = m(V’(t) ,x(t) ) = sup H(^,x,u) . 

uGU 

(b) At the terminal time tj , the relations 

(27) < 0, m(V’(t 1 ) ,x(tj) ) = 0 



are satisfied. It turns out that if V’(t) , x(t) and u(t) satisfy 
systems 



(28) 



dx i _ <9H 
dt dxj 



i = 0,1,. . . ,n , 



(29) 



dt 




0 , 1 , . . . , n 






and condition (a), then the time functions tp 0 (t) and m ( rfi (t ) , x (t ) ) 
constant. Thus (27) may be verified at any time t, t Q < t < tj , 
not j ust t x . 



to 



are 

and 
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F o rmu 1 at ion 



The Blaschke-Lebesgue theorem states that the Reuleaux triangle has 
the least area of all plane convex sets K of the same const nat width 

b. The minimum area is — . In the Following, without loss oF 

generality, we assume b = 1 . We also write our Formulas concerning a 
plane convex body K without speciFy ing K. For example (6) will be 
written as 

(30) h(0) + h(0) = p(0) . 

Thus we can rewrite Formulas (5), (0), (7), (10) and (11) as Follows: 



(31) 


A = 


1 C 2 * 

k h (0)p(0) 60 

J o 


(32) 


h(0) 


+ h(0) = P(6) , 


(33) 


h(0) 


+ h(^+^) = 1 , 


(34) 


p(0) 


+ p(0+ir) = 1 , 


(35) 


0 < 


P(0) < i • 



We 

can be 



also have the Following Formula For any convex body K which 
derived From Frenet-Se r ret Formulas (Flanders [12]). 



h(0) = x(0) • t(0) , 



where x denotes 
point where the 
Through the 
width there are 



the position vector and t unit tangent vector at the 
outward unit normal is given by (cos 0, sin 0). 
endpoints oF any diameter oF a set K oF constant 
support lines oF K perpendicular to that diameter. 
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Furthermore there is a diameter for which the corresponding support 
lines are tangent to the curve (Eggleston [11] , P. 126). Let one 
endpoint of such a diameter be taken as origin. Using (33) and h(0) 
we obtain 



(36) 


h (0) = 1, h ( 7r ) = 0 , 


(37) 


h (0) = 0 , h(jr) = 0 . 



Substitute (33) and (34) in (31) to derive 



(38) 

Let 


A = 1 f [1 + 2h(0)p(0) - h(0) - p(0)l 60 . 
J 0 


(39) 


Xl (tf) = 2h(0) - 1 , 


(40) 


x 2 (0) = Xj(0) = 2h(0) , and 


(41) 


u(6») = 2 p(0) - 1 . 



If we now substitute (39), (40) and (41) in (32), (35), (36), (37) 

and (38) , we obtain the following approach to proving the Blaschke- 
Lebesgue theorem: 



(42) 

subject to 


f7T 

Minimize 1 (1 + >: 1 (0)u(0)) d0 , 

J 0 


(43) 


*1 = x 2 , 


(44) 


X 2 = u - X x , 


(45) 


X 

o 

'w' 

II 

X 

h-* 

/— S 

II 
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(46) x 2 ( 0 ) = 0^ x 2( 7r ) = 0 ' 

(47) I u | < 1 . 



Anal ys i s 



The B 1 aschke selection theorem states that every infinite 
sequence of closed convex subsets of a bounded portion of R n contains 
an infinite subsequence that converges to a closed nonempty subset 
(Eggleston [11] 9 P. 64) . The Blaschke selection theorem impl i es that 
the minimum area exists. 

Minimizing (42) subject to (43)-(47) is equivalent to minimizing 

(48) [ x 1 (fl)u(fl) dO 

J 0 



under the same constraints. We proceed to use Pontryagin ’s maximum 
principle to minimize (48) . 

Let 



(49) 



dx Q 

66 



Xi(0)u(0) . 



Let, t/>(6) = ( (#) , V’i (6) , rp 2 (6) ) be the auxiliary vector. Use (25) to 

obta i n 



(50) 



H = ^q x i u + ^i x 2 + V> 2 C u — x l ) = ( V , o x i + V’ 2 ) u + ^i x 2 “ t' ! 2 x i 



We then use (24) to write down differential equations for v 



(51) 



d^o _ 

d<? ” 



AH 

dx Q 



- 0 



(52) 



d tpi 
66 ~ 



ffi = 



11 



(53) 



dt/> 2 

d0 



d\\ _ 

5x 2 



-V>i 



There exists el nonzero continuous vector 0(0) satisfying (51)-(53) 
such that 

(54) max H (x (0) , u , 0(0) ) = H (x (0) , u (0) , 0(0) ) . 

M<1 



Now we consider cases to analyze the maximum of H given in (50) 
as a linear function of u. 

Case 1 (a.) : rp 0 ^ 0 and ^ 0 x i + V’2 ¥■ 0 (not identically equal to 

zero) lead to a cont rad i ct i on . 

In this case 

f 1 if ^0*1 + > 0 9 

(55) u = l 

l-l if V’o x i + ^2 < 0 . 

Using (43) and (44) we obtain x 1 = u -f A sin (0 -q) for some A and a. 
Continuity of V’o ■* x i a ^d \l } 2 imply that V’o x i + ^2 — ( u + A sin (r-a)) 

+ V ; 2 * s continuous at a switching point r. Continuity at r implies 

(56) 4- jpQ A sin (r-a) + V’ 2 ( r ) = “V’o + *l>o A sin(r-a) + ^ 2 ( r ) 

Hence using (56) we conclude that rp 0 = 0 which is a contradiction. 
Case 1(b) : Assume ip Q £ 0 and V>o x i + ^2 = 0. 

In this case ^2 = ”V , o x i • Hence — - V’o x i = — _ ^o x 2 ? an d 

(57) - V’ou = V’o x 2 = V , o( u " x i) • 

Hence 



-V’oXi - V’ou = V’ou - t/’o x i > 
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which implies u = 0. Using (41) we 
corresponds to a circle of radius ~ 
than the minimum. 



Case 2 . 


Suppose rp 0 


= 0 


• 


In this 


case 


^2 £ 


0 . 


Since i/> 2 


we know that 


V’ is 


a nonzero 


vector . 




r 


i 


if 


nT 

V 

C 


(58) 


u = { 


-i 


if 


N) 

A 

C 


D i f f e rent i al 


equat i ons 


(52) 


and (53) 


(59) 


d6» “ 


i’ 2 ' 






(60) 


II 

lb 


-V’i . 







conclude p(0) = — which 

giving the maximum area rather 

= 0 and (53) imply = 0 . But 
Hence using (55) we obtain 



reduce to 



An analysis similar to Pontryagin [19, pp . 2 / -35] will give the 

switching curve in Figure 2 where 



(61) 




above the cuve 
below or on the curve 



We can now interpret the system of differential equations (43) 
and (44) as equations of motion. Our objective would be to get from 
(1,0) to (-1,0) in such a way that we minimize (48). One can show 
that this is equivalent to minimizing the action of the moving 
ob j ect . 

Optimal trajectory starts at (1 ,0) . Since (1 ,0) is above the 
curve in Figure 2, we use u = -1 until the trajectory intersects the 
switching curve on the x 2 axis. We will then use u = 4-1 which will 
lead to (-1,0). If we interpret 0 as time, then the total time from 
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(1,0) to (-1,0) is We realize that x 1 = 1 and x x = -1 are 

solutions o-f (43)- (46) in neighborhoods of 6 = 0 and 6 = n 
respectively. Hence in order to use the total time suggested by 
functional (48), the object starts at (1,0) and waits there for time 
a < Then the object goes to (0,->|~3) and then to (-1,0). During 

the time interval (^-fa,7r) the object waits at (-1,0) for a period 
°f §-a. In this way the total time would be 7r . See Figure 3. 

We use differential equations (43)-(46) and (62) below to obtain 
(63) . 



(62) 



(63) 



u 



1 , 0 < 0 < a 

- 1 , a < 0 < Z+a 

l 

1 , ^ + a < 0 < ^ + a 

-1 , a + ^ < 0 < ir . 

I J 



h (6> 



1 + Xj(0) 

2 



1 , 






cos (0-a) , 

1 + cos (0 + ^ - a) , 

0, 



0 < 0 < a 
a < 0 < ^ + a 
g + a < 0 < ^ + a 
o + ^ < 0 < ir . 



However (63) gives the support function of a Reuleaux triangle of 
width 1. The angle a corresponds to the fact that a rotation of a 
Reuleaux triangle by angle a will result in a new support function 
where 6 is replaced by 0 < a < g. We can now calculate (42) to 

obtai n 
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(64) 



Minimum area = -1 

4 



(l + xjC^uCfl)) do = 1 - r 

0 



Hence by comparison of case 1 and case 2 we choose the minimum, which 

7T - -[3 



1 s 



as desired 
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Figure 1. Reuleaux Triangle 
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+ 1 on or below the curve 
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Figure 2. Switching Curve 
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